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Using numerical simulations we examine colloids with a long-range Coulomb interaction confined in a two-dimensional trough 
potential undergoing dynamical compression. As the depth of the confining well is increased, the colloids move via elastic 
distortions interspersed with intermittent bursts or avalanches of plastic motion. In these avalanches, the colloids rearrange 
to minimize their colloid-colloid repulsive interaction energy by adopting an average lattice constant that is isotropic despite 
the anisotropic nature of the compression. The avalanches take the form of shear banding events that decrease or increase the 
structural order of the system. At larger compressions, the avalanches are associated with a reduction of the number of rows 
of colloids that fit within the confining potential, and between avalanches the colloids can exhibit partially crystalline or even 
smectic ordering. The colloid velocity distributions during the avalanches have a non-Gaussian form with power law tails and 
exponents that are consistent with those found for the velocity distributions of gliding dislocations. We observe similar behavior 
when we subsequently decompress the system, and find a partially hysteretic response reflecting the irreversibility of the plastic 
events. 


1 Introduction 

Collectively interacting colloidal particles are often used as 
model systems to investigate various features of equilibrium 
and non-equilibrium phenomena-. Due to their size scale, 
colloids provide the advantage that microscopic information 
on the individual particle level can be directly accessed, some¬ 
thing which is normally difficult or impossible in smaller scale 
systems such as nanoparticles, molecules, or atoms^. Ad¬ 
ditionally, there are various methods such as optical tech¬ 
niques^ for controlling colloidal ordering and manipulating 
individual colloids. Examples of phenomena that have been 
studied with colloids include two-dimensional melting transi¬ 
tions^, solid-to-solid phase transitions^, glassy dynamics^, 
commensurate and incommensurate phases^ - —, depinning be¬ 
haviors-^ - —, self-assemblyiii^, and dynamic sorting— - - 1 2 * . It 
is also possible to use colloids to study plastic deformation 
under shear in crystalline^ or amorphous^ colloidal assem¬ 
blies. In crystalline materials, plastic deformations occur via 
the motion of dislocations which often occur in bursts of ac¬ 
tivity or avalanches^"- 8 . Certain studies that may be difficult 
to undertake in other systems become feasible to perform with 
colloidal systems, such as observations of changes in the par¬ 
ticle configurations and dynamics during compression. Due to 

1 Theoretical Division, Los Alamos National Laboratory, Los Alamos, New 
Mexico 87545, USA. Fax: 1 505 606 0917; Tel: 1 505 665 1134; E-mail: 
cjrx@lanl.gov 

2 Department of Physics, Wabash College, Crawfordsville, Indiana 46556, 

USA. 


the relative softness of charge-stabilized colloidal assemblies, 
compressional studies can be performed in this system over a 
wide range of parameters and packing densities. In recent ex¬ 
periments, Varshney et al . 29 demonstrated that it is possible to 
create a quasi-two-dimensional colloidal raft system confined 
between two barriers, and to then compress or decompress 
the raft in order to dynamically change the packing density. 
In these experiments the colloid-colloid interaction was more 
complex than a simple repulsion, so that transitions from a 
loose-packed amorphous solid to a much denser but still amor¬ 
phous solid were observed; however, during the compression, 
various plastic rearrangements occurred that produced hys¬ 
teresis across the compression and decompression cycled. 

Several studies have addressed order-disorder transitions of 
repulsive colloids confined between two parallel walls, where 
the transitions are associated with changes in the number 
of rows of colloids that fit between the walls. Crystalline 
states arise for integer numbers of rows, while disordered 
states occur when partial rows, or a non-integer number of 
rows, are present 30,31 . It should be possible to create a two- 
dimensional (2D) system of repulsively interacting colloids 
that tend to form a triangular lattice, confine the system by 
barriers or an anisotropic trapping potential that can be dy¬ 
namically changed to compress the colloidal assembly in one 
direction, and observe how the configuration changes under 
compression as the colloids rearrange in order to minimize the 
repulsive colloid-colloid interactions. For example, in numer¬ 
ical studies of confined charged particles in anisotropic traps, 
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the particle configuration changes as a function of the trap 
width— - —. In experiments on ion trapping systems where the 
ions are confined in a quasi-one-dimensional potential with a 
tunable strength, structural changes in the ion configuration 
were observed as the assembly was compressed and decom¬ 
pressed, such as a transition from a single row of ions to a 
zig-zag or kink stat e 35 ^ 6 . It should also be possible to create 
dynamical confining potentials or barriers for dusty plasma 
crystals 37 . 

In this work we numerically study a 2D system of col¬ 
loids with long-range Coulomb repulsive interactions placed 
in a quasi-one-dimensional confining potential such that as 
the depth of the potential is increased, the colloidal assembly 
is compressed in one direction. Since the colloids can mini¬ 
mize their interaction energy by adopting a triangular ordering 
with equal lattice spacing in all directions, the colloids adjust 
their positions to make the local lattice structure as isotropic 
as possible during compression. Colloidal motion occurs both 
through slow elastic distortions and via abrupt avalanches in 
which plastic rearrangements occur. During the avalanche 
events, the average distance between nearest neighbor colloids 
increases. Dislocations can be created or annihilated at the 
sample surface during the plastic events, causing the colloids 
to move along shear bands. For high compressions, we find 
that avalanches are associated with reductions of the number 
of rows of colloids that fit in the confining potential. Between 
the row reduction events, the colloids can adopt partial crys¬ 
talline or smectic order. We also find that the colloid velocity 
distribution during the avalanche events is non-Gaussian with 
a power-law tail, and that the exponents are consistent with 
those observed for velocity distributions in 2D dislocation sys¬ 
tem s 32 i 38 i 39 . Under decompression, we observe a similar set 
of dynamics and find some hysteresis between the compres¬ 
sion and decompression cycles; however, there are no large- 
scale hysteretic effects since the particle-particle interactions 
are purely repulsive. 

2 Simulation and System 

We consider a two-dimensional assembly of colloidal parti¬ 
cles interacting repulsively via a long-range Coulomb poten¬ 
tial. We employ periodic boundary conditions in both the v 
and y-directions for a system of size L x xL y . The sample con¬ 
tains N = 256 colloids and a single trough potential which 
produces a force Ff = F p cos(2nxi/L x )x. As F p is increased, 
the particles are forced closer together in the v-direction. The 
dynamics of a single colloid i are obtained by integrating the 
following overdamped equation of motion: 

JT* . N 
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Fig. 1 Under compression, achieved by increasing F p , the particle 
motion occurs through both slow elastic distortions as well sudden 
plastic events or avalanches that are associated with spikes or jumps 
in the following quantities: (a) N~ l dE/dF p , the change in the 
particle configuration energy, vs F p . (b) (d m ), the average distance 
between nearest-neighbor particles, vs F p . Inset: a blowup of the 
main panel showing that plastic events are associated with increases 
in {d m ). (c) F x (blue) and F y (green), the average transient forces on 
the particles in the x and y directions, vs F p . (d) P { g, the fraction of 
six-fold coordinated particles, vs F p . Grey bands across all panels 
correspond to the intervals of F p from which the particle trajectories 
in Fig. |2] are taken. 
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Here T] is the damping constant, R^-) is the position of par¬ 
ticle i(j), Rij = \Rj — R;|, and the particle-particle interac¬ 
tion potential is V(Rij) = q 2 Eo / R ^, where Eo = Z* 2 /47l££oao, 
q is the dimensionless interaction strength, Z* is the effec¬ 
tive charge of the colloid, and £ is the solvent dielectric con¬ 
stant. We treat the long-range image interactions using a real- 
space Lekner summation method—. Lengths are measured in 
units of ao, time in units of T = 77 / Eq, and forces in units 
of Fq = Eo/ao. We increase F p in small increments of 0.001 
over the range F p = 0 to 10 in the first part of the work, and 
in increments of 0.01 from F p = 0 to 100 in the second part 
of the work. During each increment we measure the average 
transient force response on each particle in the two directions 
Fx = if lf4>{*(')•* and F y = if F^r) -y, along 
with the changes in the particle-particle interaction energy 
E = Y 4 V(Rij)* the total energy of the system E tot = E + 

Lf F p sm(2nxi/L x ), local ordering P 6 = AT 1 5 ^ S(zt - 6), 
where Zi is the coordination number of particle i obtained from 
a Voronoi tesselation, and (d m ) = N~ l d l m , where d l m is the 
spacing between particle i and its nearest neighbor as identi¬ 
fied from a Voronoi tesselation. 

3 Compression 

In Fig. QJa) we plot the change in the energy of the parti¬ 
cle configuration N~ l dE/dE p versus F p . Figure [ljb) shows 
the corresponding (d m ) 9 the average spacing between nearest- 
neighbor particles, Fig. |TJc) shows the cumulative transient 
forces F x and F y , and Fig.HJd) shows the fraction of six-fold 
coordinated particles P { 5. In Fig. [2] we plot the particle posi¬ 
tions and trajectories over intervals of 8 F p = 0.1 which corre¬ 
spond to the gray shaded bands in Fig.Q] We find two distinct 
types of behaviors. Elastic distortions occur in the smoothly 
changing portions of dE/dE p and (d m ), in the same regions 
where F x and F y are close to zero. Sudden plastic events or 
avalanches are associated with peaks in dE/dE p , jumps in 
(d m ), and spikes in F x and F y . The particle configuration en¬ 
ergy E increases with increasing E p as the repulsively interact¬ 
ing colloids are forced together, which accounts for the overall 
positive value of dE/dF p . The cusp in dE/dF p near E p = 1.0 
appears when F p becomes large enough that the particles can 
no longer spread to cover the entire sample but instead segre¬ 
gate into the bottom of the potential, with a particle-free region 
appearing at the potential maximum that becomes wider with 
increasing F p . 

The spikes in dE/dE p generally have negative values, indi¬ 
cating that they are associated with drops in the configuration 
energy. The plastic events are also associated with an increase 
in the average spacing between nearest-neighbor particles, as 
indicated in Fig. |TJa,b) where negative spikes in dE/dE p cor¬ 
relate with positive jumps in (d m ). The inset in Fig.HJb) shows 



Fig. 2 The particle positions (colored dots) in the potential as the 
particle assembly is compressed along the x direction by increasing 
F p . Blue particles are stationary or slowly moving while red 
particles are moving the most rapidly. Small black dots indicate the 
particle trajectories over a force interval of 8F p = 0.1. (i) At 
F p = 2.0, the motion is predominantly elastic, (ii) At F p = 2.5, an 
avalanche occurs in the form of shear bands, (c) At Fj = 4.8 the 
avalanche activity is more localized, (d) At F p = 6.6 only elastic 
motion occurs, (e) At F p = 7.2 there is a large avalanche, (f) At 
F p = 8.0 the motion is predominantly elastic. 
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Fig. 3 (a,d) Particle positions colored by the amount the particle moved during an increment of F p , with blue particles nearly stationary and 
red particles moving the largest amount. (b,e) Velocity distributions P(V X ), P(V y ), and P(|V|), from top to bottom. The same color scale 
shown for P(|Vj) is used in panels (a) and (d). (c,f) P(|V X |), P(|V^|), and P(|V|) plotted on a log-log scale. The solid line in panel (c) is a 
power-law fit P(|V|) oc \V\ T with an exponent of T = —2.2. Panels (a,b,c) are for a force interval of 3.4 < F p < 3.5 where an avalanche occurs, 
and panels (d,e,f) are for a force interval of 3.6 < F p <3.7 where the behavior is elastic and the distributions are much sharper. 


a blowup of the plastic events near F p = 7.2 to better highlight 
the discrete nature of the jumps in (d m ). Since the colloid- 
colloid interactions are repulsive, the configurational energy 
increases under compression as the particles elastically ap¬ 
proach each other along the compressed direction, while dur¬ 
ing the plastic avalanches, the closest spacing between parti¬ 
cles increases so that the configurational energy is reduced. 
In Fig. EC) the peaks in F x coincide with peaks in F y , in¬ 
dicating that the particle motion is occurring almost equally 
in both the x and y directions during an avalanche. Fig. [ljd) 
shows that initially P^ ^0.95 at F p = 0.0 where the system 
forms a slightly distorted ordered triangular lattice. As F p in¬ 
creases, there are numerous very small peaks in P§ caused by 
the Voronoi construction having difficulties with the sample 
edges; however, there are several real larger scale jumps in P& 
that are correlated with the avalanche events. The increases 
and decreases in Pg occur when dislocations, which have five¬ 
fold or sevenfold rather than sixfold coordination, are created 
or annihilated during an avalanche. 

In the avalanche illustrations of Fig. [2 red particles move 
significantly during the compression increment while blue 
particles experience little to no motion. Figure [2a) shows 
the state near F p = 2.0 for an interval in which there are no 
spikes or jumps in the quantities plotted in Fig. Q] indicating 
that the particles are undergoing elastic motion. Here the parti¬ 
cles are arranged in twelve vertical rows. The avalanche near 


F p = 2.5 is shown in Fig. (2b) where the particles transition 
from twelve to eleven vertical rows. Here, motion occurs in 
localized bands, indicative of a shear banding effect. These 
bands generally form zig-zag patterns, and there are also re¬ 
gions undergoing local rotation. Figure[2c) illustrates another 
avalanche interval near F p = 4.8, where similar shear banding 
motion appears. In this case, there is considerable motion of 
the particles on the outer edges of the sample and the particles 
transition from ten to nine vertical rows. Figure [2d) shows 
another elastic region near F p = 6.6 where there is little mo¬ 
tion. Near F p = 7.2, a large plastic event occurs as illustrated 
by the bands of motion in Fig. |2e) as the system transitions 
from nine to eight rows. For 7.2 < F p < 10, the system be¬ 
haves elastically, as shown in Fig.[2f) for F p = 8.0 where lit¬ 
tle motion occurs. Further compression of the system beyond 
F p = 10 is described in Section [4] 

3.1 Velocity Distributions 

We next examine histograms of the particle velocities at dif¬ 
ferent values of F p . Figure [2a) shows the particle positions 
for 3.4 < F p < 3.5, where the color denotes the average speed 

|Vj = yjV} + V y of each particle, highlighting the heteroge¬ 
neous displacements of particles in bands. In Fig. [2b) we 
plot P(V X ), P(V y ), and P(|Vj) for the same force interval. 
The distributions are obtained by averaging the velocities of 
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Fig. 4 A heightfield plot of P(|V |) vs F p , where as indicated by the 
scale bar blue values of \V | occur with high frequency while red 
values of \V | occur with low frequency. The total height of each 
vertical column indicates the maximum range of \V | for that value of 
F p . The maximum value of this range gradually increases with 
increasing F p . 
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Fig. 5 Fxy Fy N dE jdF p , and P ^ for a system subjected to a 
maximum compression of F p = 100. The density of avalanches is 
much lower for F p > 10 than in the F p < 10 regime already 
discussed, and the avalanches at higher F p are associated with a 
partial or complete reduction in the number of rows of particles in 
the system. 


each particle over ten subintervals of 8 F p = 0.01. The plas¬ 
tic motion that occurs during the avalanche event produces 
non-Gaussian velocity distributions, as shown more clearly in 
Fig. Etc) where we plot P(|V*|), P(|Vy|), and P(|V|) on a log- 
log scale. The solid line is a power law fit of P(\V\) oc |y| T , 
where z — —2.2. In Fig. EId,e,f) W e show the same quanti¬ 
ties for 3.6 < F p < 3.7 which corresponds to an interval in 
which the response is elastic, as indicated in Fig. Etd). The 
width of the velocity distributions in Fig. Ete) is significantly 
smaller than for the plastic flow event illustrated in Fig. Ob). 
Due to the sharpness of the distributions, the plot of P(\V |) on 
a log-log scale in Fig. Of) cannot be fit to a power law. In 
Fig-Og>h,i) we show the same quantities for 7.1 < F p < 7.2, 
corresponding to an interval containing an avalanche as high¬ 
lighted in Fig. [T] Here the motion follows a banding pat¬ 
tern. The velocity distributions have strongly non-Gaussian 
features, as shown in Fig. Ogd 1 ). A fit of P(|V|) oc |y| T in 
Fig. Oi) gives T = —1.9. For 7.3 < F p < 7.8 there are few 
plastic rearrangements and the velocity distributions are very 
similar to those shown in Fig. Oc,f) for the F p = 7.0 case. 
In general, we find that during the avalanches, the velocity 
distributions of the particles have a non-Gaussian form with 
a power law tail that can be fit with 1.9 < T < 2.6, while in 
the elastic regimes, the width of the velocity distributions is 
strongly reduced. 

Previous studies of sheared crystalline materials showed 
that avalanches are associated with the motion of disloca¬ 


tions^ " — , and that various quantities such as the dis¬ 
location velocities are power-law distributed in the avalanche 
regime. We consider the particle velocities instead of the dis¬ 
location velocities in our system; however, due to the partial 
crystalline order of our sample, avalanches are generally as¬ 
sociated with the motion of dislocations. In the dislocation 
dynamics studies, the high velocity tails of the velocity distri¬ 
butions can be fit with power law exponents of z = -2.5— and 
3.(>™ Recent computational and theoretical work showed 
that in 2D, a single dislocation has a power-law distributed 
velocity with z = — 2, while when collective effects are im¬ 
portant, larger exponents of z = —2.4 appear—. In our case, 
as dislocations move through the system, individual particles 
temporarily translate along with the dislocations, so that over 
sufficiently short time windows the overall form of the veloc¬ 
ity distributions of the particles and the dislocations should be 
similar. It is difficult to extract the exact exponents for our sys¬ 
tem due to the finite width of the sample and the fact that some 
avalanches are associated with only one moving dislocation 
while others contain multiple moving and interacting disloca¬ 
tions. The width of the distribution is also dependent on the 
packing fraction. As the particle density increases, the maxi¬ 
mum possible particle velocity also increases. This is shown 
in Fig. 0] where we plot a color height map of P(\V |) versus F p 
for intervals of 8 F p = 0.1, highlighting the maximum extent of 
the high-velocity tails of the distributions. As F p increases, the 
particle packing density increases and the maximum achiev¬ 
able particle velocity also increases significantly. 
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Fig. 6 (a) Real space image of particle locations in a portion of the 
sample at F p = 10 showing that there are eight rows of particles, (b) 
The corresponding structure factor S( k) at F p = 10 indicates strong 
triangular ordering, (c) Particle positions and trajectories in a 
portion of the sample during an avalanche event at F p = 11.6. (c) 
S(k) during the avalanche shows a disordered structure, (e) Real 
space image of the particles in a portion of the sample at F p = 11.8 
where there are now seven rows of particles, (f) The corresponding 
S(k) shows triangular ordering. 







(c) 



Fig. 7 (a) Real space image of particle locations in a portion of the 
sample at F p = 39, where there are six rows of particles, (b) The 
corresponding S(k). (c) Real space image of particle locations in a 
portion of the sample after an avalanche at F p = 40. There are now 
five rows of particles, (d) The corresponding S( k) is more smeared, 
indicating a smectic-like structure. 
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Fig. 8 (a) Real space image of particle locations in a portion of the 
sample at F p = 80, where there are five rows of particles, (b) The 
corresponding S( k). (c) Real space image of particle locations in a 
portion of the sample after an avalanche at F p = 87. There are now 
five rows of particles, (d) The corresponding S(k) has a smectic 
structure. 


In Fig. [5] we plot simultaneously F x , F y , the derivative of the 
total energy dE ioi /dF p , and Ffc versus F p as the sample is com¬ 
pressed all the way to F p = 100. The density of avalanche- 
induced peaks is much smaller for F p > 10 than for F p < 10, 
while the avalanche jumps at higher F p values clearly coin¬ 
cide with sharp changes in P { 5 . In general, plastic avalanche 
events that occur for F p > 10 are associated with a partial or 
complete reduction of the number of rows of particles that fit 
across the potential well. In Fig. 0a) we show the particle 
positions in a portion of the sample at F p = 10, with the cor¬ 
responding structure factor S( k) = N~ l | e~ lk ' Ti | 2 plotted in 
Fig. 0b). Here, eight rows of particles fit across the poten¬ 
tial well and the system has significant triangular ordering as 
indicated by the six-fold peaks in S(k). Figure 0e,f) shows 
the real space particle positions and S(k) at F p = 11.8, where 
there are only seven rows of particles that still have mostly tri¬ 
angular ordering. Figure 0 shows that between F p = 10 and 
F p = 11.8 a plastic avalanche event occurs, as indicated by 
the spikes in F x , F y , and dE ioi /dF p . We illustrate the corre¬ 
lated bands of particle motion that occur during the F p = 11.56 
event in Fig. 0c). Figure0d) shows that during the avalanche, 
the sixfold peaks in S( k) are heavily smeared, indicating the 
disordering of the structure as dislocations move through the 
system. We observe similar S( k) signatures for avalanches at 
higher F p . In some cases, the number of rows is not reduced 
uniformly through the entire system in a single avalanche; in¬ 
stead, a portion of the system collapses to the smaller num¬ 
ber of rows in one avalanche event, followed by a second 
avalanche event at a slightly higher value of F p that completes 
the row reduction process. 

In Fig. 0a, b) we show the real space and S( k) plots for 
a sample at F p = 39 where there are six rows of particles. 
Here there are patches in the system where the triangular or¬ 
dering has been distorted into a rhomboidal ordering, as seen 
by the nonuniform spacing between the peaks in S(k). Fig¬ 
ure 0 shows that an avalanche occurs near F p = 40. After the 
avalanche, at F p = 44, the real space plot in Fig. 0c) has only 
five rows of particles, indicating that the avalanche was as¬ 
sociated with a row reduction. The structure in Fig. 0c) has 
a smectic characteristic as indicated by the one-dimensional 
smearing of the peaks along the k x direction in the S( k) plot 
in Fig.0d). For 42 < F p < 80, the system retains five rows of 
particles; however, there are still avalanche events associated 
with the system shifting between different partially crystalline 
structures that still contain five rows. For example, Fig. 0a) 
shows a rhomboidal ordering of the five rows of particles at 
F p = 80, and Fig. 0b) shows the corresponding S(k). In the 
range 80 < F p < 87, several avalanche events occur that are 
associated with a transition from five to four rows of parti¬ 
cles. In Fig. 0c, d) we show the particle positions and S( k) at 
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Fig. 10 (a,d) Particle positions colored by the amount the particle moved during an increment of F p , with blue particles nearly stationary and 
red particles moving the largest amount. (b,e) Velocity distributions P(V X ), P(V y ), and P(|V|), from top to bottom. The same color scale 
shown for P(|V|) is used in panels (a) and (d). (c) P(|\^|), P(|\^,|), and P(|V|) plotted on a log-log scale. Solid line is a power-law fit of 
P(\V\) oc \V\~ T , where T = —2.0. Panels (a,b,c) are for a force interval of 5.9 > F p > 5.8 during decompression, where a large plastic event 
occurs. Panels (d,e) are for a force interval of 5.8 > F p > 5.7 during decompression, where the system behaves elastically. 
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Fig. 9 Hysteretic response measurements. Blue curves are obtained 
during the initial increasing sweep of F p , and green curves are from 
the decreasing sweep in which F p is reduced back to zero. Grey 
bands are at the same locations as in Fig.Q] and indicate the 
intervals during the increasing sweep of F p at which the particle 
trajectories in Fig. [2| were obtained, (a) P { 5 vs F p . The blue curve for 
increasing F p also appears in Fig.QJd). (b) F x vs F p . The blue curve 
for increasing F p also appears in Fig. 0c). ( c ) ®xx vs F p . 


F p = 87, where there are four rows of particles in a smectic 
arrangement. 

5 Decompression and Hysteresis 

To explore hysteretic effects, once the system is fully com¬ 
pressed to F p = 10 we gradually decompress the sample by 
decrementing F p back to zero, using the same sweep rate for 
F p as during the original compression. We find that, just 
as during compression, under decompression the system un¬ 
dergoes combinations of slow elastic distortions interspersed 
with sudden avalanche rearrangement events. In Fig. 0a, b) 
we plot Pa and F x versus F p for the compression and decom¬ 
pression cycles in a sample where the maximum compres¬ 
sion value is F p = 10. Figure 0a) shows that the Pa curves 
for compression and decompression overlap completely above 
F p ~ 8.5, while in Fig. 0b) there are no avalanche events in 
this same force interval, indicating that the system is behav¬ 
ing reversibly. Avalanche events are completely suppressed 
during the decompression until F p has decreased to F p = 6.0, 
where the first large peak in F x appears in Fig. 0b) accom¬ 
panied by a large dip in Pa in Fig. 0a). During this event, 
a number of dislocations enter the system and the average 
spacing between the colloids increases. In general we ob¬ 
serve fewer plastic events during the decompression cycle than 
during compression; however, the plastic events that do occur 
during decompression tend to be larger. Two large plastic de¬ 
compression events appear at F p = 4.0 and F p = 2.6, marked 
by spikes in F x and substantial drops in Pa. We find no sys¬ 
tematic trend in the hysteresis such as Pa or F x always being 
higher for one direction of the compression cycle. Instead, 
the plastic avalanches occur at different values of F p for each 
driving direction, indicating the irreversibility associated with 
the plasticity in this system. The lack of systematic hystere¬ 
sis is different from the behavior observed for compression- 
decompression cycles in colloidal raft experiments^. In those 
experiments, the particle-particle interactions had both repul¬ 
sive and attractive components, so that during the initial com¬ 
pression the particles could be captured by the attractive por¬ 
tion of the interparticle potential. In contrast, in our system all 
of the particle-particle interactions are smoothly repulsive. 

Figure 0c) shows a**, an element of the stress tensor: 

G xx = ^2 (Fij)x(?ij)x ( 2 ) 

^ i j<i 

plotted versus F p for both compression and decompression. 
Here F^ is the interparticle force between particles i and j and 
fij is their relative separation. There is no significant differ¬ 
ence between the compression and decompression values of 
G xx over the range F p > 8.0, but we find significant hysteresis 
for F p < 8.0. For F p < 1.8, the < 7 ^ signal is noisy under both 
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Fig. 11 Images illustrating the hysteretic behavior of the sample 
shown in Fig. |9jc) over the range F p = 1.1 to F p = 1.8 where the 
system transitions from diagonal to vertical rows of particles. 
Particle positions are colored by the amount the particle moved 
during an increment of F p , with blue particles nearly stationary and 
red particles moving the largest amount. Small black dots indicate 
the particle trajectories over a force interval of 8F p =0.1. In (a-c) 
the system is being compressed while in (d-f) the system is being 
decompressed. (a,d) At F p = 1.1 the rows of particles are not 
oriented along the y axis. (b,e) At F p = 1.6 the rows in the 
compressing system (b) are vertical while those in the 
decompressing system (e) are diagonal. (c,f) At F p = 1.8 the rows 
are vertical during both compression and decompression. 


compression and decompression. At very low F p , the parti¬ 
cles are spread out across the entire sample, with no gap. At 
F p = 1.5, a particle-free gap opens at the location of the po¬ 
tential maximum, and the particles form rows that are rotated 
by a finite angle from the y direction, as shown in Fig. fTlT a.f). 
This diagonal ordering reduces the energy of the system by 
minimizing the compression of the lattice, but when F p > 1.8 
this small energy advantage is lost and the particles transi¬ 
tion as illustrated in Fig. M into a vertical row structure of 
the type shown in Fig.QjJc) that is aligned with the direction 
of the confining potential. For 1.8 < F p < 8.0, q** remains 
hysteretic but is no longer noisy; instead, the curve is smooth 
with distinct jumps, indicating the hysteresis in the value of F p 
at which the number of rows in the sample is reduced during 
compression or increased during expansion. The row tran¬ 
sitions are marked by labeled arrows in Fig. 0c), where we 
observe an average offset of A F p ~ 0.7 between the reduction 
in the number of rows from n to n — 1 upon compression and 
the increase in the number of rows from n — 1 to n upon de¬ 
compression, with the transition during compression falling at 
a higher value of F p . For example, the 11-10 transition dur¬ 
ing compression occurs at F p = 3.2 while the 10-11 transition 
during decompression occurs at F p = 2.7. 


During decompression, the particle velocity distributions in 
the elastic regimes and during avalanches have the same fea¬ 
tures described above for compression. Figure[l3a,b,c) shows 
the particle configuration, velocity distributions, and log-scale 
plot of P(|V X |), P(\V y \), and F(|V|) during the decompression 
cycle over the range 5.9 > F p > 5.8. The velocity distribu¬ 
tions are clearly non-Gaussian, and the power-law fit of P(\V |) 
in Fig. flQl c) has an exponent of T = —2.4. Figure ITQF d.e.f) 
shows the same quantities over the range 5.8 > F p > 5.7 where 
elastic motion occurs. Here the velocity distributions are much 
sharper. 

6 Summary 

We have numerically investigated a monodisperse assembly of 
repulsively interacting Yukawa colloids undergoing dynami¬ 
cal compression in a trough potential. During the compres¬ 
sion the motion is predominately elastic as the particles grad¬ 
ually move closer together. This elastic motion is interspersed 
with occasional sudden plastic rearrangements or avalanches 
in which the colloids shift position in order to increase the av¬ 
erage spacing between neighboring particles. The avalanches 
take the form of local shear banding events in which dis¬ 
locations can be created or annihilated. During the plastic 
events, the colloidal velocity distributions are non-Gaussian 
with power law tails that have exponents ranging from T = 
—1.9 to T = — 2.5. This is consistent with the velocity distribu¬ 
tions found for moving dislocations during avalanche events. 
For larger compressions, the avalanches are generally associ¬ 
ated with row reduction events where the number of rows of 
colloids that fit inside the trough is partially or completely re¬ 
duced by one. After these avalanche events, the colloids can 
exhibit partial crystalline or smectic type structures. During 
decompression, we observe similar avalanche behaviors; how¬ 
ever, the avalanche events do not occur at the same values of 
substrate strength, indicating the occurrence of irreversible be¬ 
havior. This system could be experimentally realized by con¬ 
fining colloids or other charged particles such as dusty plas¬ 
mas in anisotropic traps where the barriers or the trap width 
can be changed dynamically. 
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